International Journal of Power Electronics and Drive System (IJPEDS) 
Vol. 11, No. 1, March 2020, pp. 220~227 
ISSN: 2088-8694, DOI: 10.11591/ijpeds.v11.11.pp220-227 g 220 


Structural correction of nonlinear dynamics of frequency- 
controlled induction motor drives 


Vladimir L. Kodkin, Aleksandr S. Anikin, Aleksandr A. Baldenkov 
Departement of theoretical fundamentals of electrical engineering, South Ural State University, Russia 


Article Info 
Article history: 


Received Jul 7,2019 
Revised Aug 12, 2019 
Accepted Nov 18, 2019 


Keywords: 


Dynamic positive feedback 
Frequency analysis 
Frequency response 
Induction motor drive 
Nonlinear dynamic 


ABSTRACT 


Frequency analysis of processes in the frequency-controlled induction motor 
drive s, as in automatic control systems with dynamic non-linearities, are 
presented in this paper. For the first time, the dynamic transfer function in 
terms of electromagnetic torque (m) from absolute slip determined by the 
load and frequency of the supply voltage are proposed. To study the dynamic 
characteristics, a method for constructing a family of frequency 
characteristics for fixed values of absolute slipss and frequencies of the 
supply voltage are described. The frequency characteristics corresponding to 
these transfer functions are constructed in the Simulink application of the 
MatLab software. Based on results of the analysis the structural correction 
methods that improve the dynamics of such nonlinear systems are proposed. 
This paper is an example of how the initial complicated (but more accurate) 
interpretation of nonlinearity allowed us to find a new best solution to the 
problem of controlling a complex dynamic object. 
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1. INTRODUCTION 

The widespread use of frequency converters for controlling asynchronous motors in recent years has 
created the impression that there are no problems in the field of induction motor drive (IMD). However, the 
results of studies [1-5] show that the dynamic characteristics of these systems require a more in-depth study 
from the point of view of the theory of control of nonlinear systems. Frequency controlled IMDs are a highly 
non-linear system. The “main” parameter determining the non-linearity of these systems is the variable 
frequency of the supply voltage. Unlike the stationary non-linearities of the regulatory systems considered in 
the 80s and 90s of the 20th century, the variable frequency in the IMD changes its frequency response. 

Variable frequency, strictly speaking, does not allow the use of a mathematical apparatus designed 
for IMDs, based on vector analysis, since the vector representation of variables over time implies the 
constancy of the frequency of the supply voltage, or the frequency of rotation of these vectors. However, due 
to the absence of another, vector methods are used in most research or educational works on IMDs [6-8], 
despite the fact that the authors quite often recognize the illegality of such an approach. 


2. PROBLEM STATEMENT 
The most effective engineering method for assessing the dynamics of electric drives is the method 
of frequency analysis [9-11]. Direct application of this method to IMDs is hampered by the presence of 
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significant non-linearities in them. The construction of the frequency characteristics of such systems involves 
a number of inevitable assumptions [12-15]. After considering the different versions of these assumptions, 
the calculation of the dynamic mechanical characteristic set forth in the Usoltsev’s monograph [12] turned 
out to be the most acceptable. The calculation is repelled by equation 1.36 on p.23 [12]. It establishes a 
connection between the current moment (m) and slip (£) at the nominal frequency @1nom: 


man t __. (1) 


= S 
(1+T2'p) [P a+r p)| +s 


Where, T, = =x — the transient time constant of the rotor, f = > — the relative slip, M;,— the critical 
2 1 
moment, Sx — the critical slip at the nominal frequency @1nom. 
At the beginning of the working characteristic (for M = 0, P > 0), the transfer function is simplified 
and reduces to a dynamic link of the 1st order: 


2M 2MrP _ 2MxR(W1-w) 


(1+T2 mE (1+T21p)Sk  (4+T2'p)Skw1 


To determine the transfer function by the electromagnetic torque (m) from absolute slip (Aw = 
w — w), we transform it to the form: 


m _ 2Męg 


WO) 5 aro (2) 


Aw (1+T2'p)Skw1 


However, the results of experiments given in [1, 16] showed that it is incorrect to extend this 
formula to all operating modes. 


3. SOLUTION 
The expression (1) allowed us to propose another variant of linearization, in which the initial 
equation takes the form: 


1 I 2M I 

m|(T,')?p? + 2Tap +1+ (£) == <B(1 + Tsp) (3) 
K k 

In this case, expression (1), the binding torque (m), relative slip (£), and quantities determined by 


the parameters of the electric motor (T’, - transition time constant; M+ sx - critical moment and critical slip, 
depending on the frequency 1) will be transformed to: 


_ 2MK(TzP+1)SKB 
~ (14+T21p)2S2 +82 (4) 


and the transfer function (2) will be represented as follows: 


2Mx(T2p+1)Sk 
w1 [(1+T21p)2 S +82] 


W(p) = (5) 


where w, — the stator supply voltage frequency. 

The working section (at 0 <m < 2M,, where M, is the rated torque) of the mechanical 
characteristics of the IMD will be described by the structural diagram shown in Figure 1. 

Equation (3) appears to be substantially nonlinear, since it will transform when the frequency of the 


supply voltage and slip of the motor change. Formula (4) can be called the dynamic Kloss formula (M = 
2'Mk 
S Sk 
Sp Ss 
p = 0, i.e. without load, formula (4) coincides with the formula given in the Usoltsev monograph [12]. These 
changing properties of the transfer function determine the change in the dynamics of IMD with frequency 
control, noted earlier in experiments. Thus, the "freezing" of the transfer functions is effective enough to 


interpret the dynamic characteristics of IMD. 





). The formulas become identical when substituted into the formula (4) p = 0, i.e. for a static state. For 
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Figure 1. Block diagram of IMD with the working section of the mechanical characteristics 


Figures 2a and 2b below show the frequency characteristics of IMD based on a low-power motor 
(0.37 kW) with a squirrel-cage rotor. To build frequency characteristics, hereinafter, the software package 
Simulink MatLab was used. The given amplitude and phase-frequency characteristics of the motor for the 
frequencies of the stator voltage of 10 Hz (Figure 2a) and 50 Hz (Figure 2b) correspond to slips at low 
(0.2M,.) - W1 and rated load (M,.) - W2. 
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Figure 2. The amplitude and phase characteristics of the IMD at a frequency of the stator supply voltage of 
10 (a)50 (b) Hz: W1 - slip at a load torque of 0.2M,., W2 - slip at a load torque of M, 


The given frequency responses well explain some problems of AC drive. Figure 2a shows that at a 
stator supply voltage frequency of 10 Hz and a small load, the phase shifts of the frequency response (W1) in 
the range from 10 to 100 rad / s vary from 40 to -90 degrees. Comparison of frequency responses at 
frequencies of stator voltage of 10 and 50 Hz shows that in the range from 10 to 100 rad / s the phase shifts of 
frequency responses have significantly different values - from 25 to -45 electric degree. This means that 
during acceleration and deceleration, the phase shifts change in such a way that a system with a stability 
margin at a frequency of 50 Hz can become unstable. The frequency of the stator voltage will be 10 Hz. This 
may be the reason for the different oscillation of the drive at different frequencies of the stator voltage, which 
were noted in [17, 18]. Thus, the nonlinearities of the transfer functions of the link of the torque generator 
Figure 1 require linearization to increase the efficiency of the electric drive and the same behavior at different 
frequencies. 

Most often, “transvector” control is used to linearize the IMD. With this control, dynamic links are 
formed in the control device that are inverse, mainly, to the inertia of the rotor circuits, with the rotor speed 
adapting to the frequency of the stator voltage, however, this does not take into account the pulsed nature of 
the frequency converter and the presence of other harmonic components in the rotor and stator currents 
(except fundamental harmonic) [19-24]. 

Another factor that significantly worsens the operation of vector control should be recognized that a 
number of parameters (for example, stator and rotor resistances) can change significantly during operation, 
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and it is almost impossible to track these changes. All this contributes to the nonlinearity of the rotor and 
stator currents during vector control and leads to a significant decrease in its efficiency [1, 4, 16]. 

Consider an alternative linearization of the transfer function Wp with respect to the electromagnetic 
torque (5). For this, in the structural diagram [25], we introduce local feedback (Figure 3 ) on the 
electromagnetic torque with the transfer function Wppp (the introduced designation will be explained below). 


2M; (Top + 1)S;, 
w [(1 + Tz'p)? SÈ + B?] 





Figure 3. Introduction of feedback Wppr 


We determine the transfer function W.gy by the electromagnetic torque: 


2M, (T2p+1)Sx 
2¢2 2 1 
_ @y[G+T2"p)*SEtBe] 2MkSk(T2p+1) 
qv j 2MkSk(T2p+1)WppF w1 [(1+T21p)Z S +62]+2MkSk(Tzp+1)WppF 
w1[(1+T21p)2SZ +62] 


We 


Under the following condition: 
wb? = —2M,S;(T2p + 1)WppFr 


That is, if the corrective element will have the following transfer function: 


=e wf? 
Wope = — reslo) (©) 


Then the transfer function Weqy by the electromagnetic torque looks like: 


W. = 2MxSx(T2p+1) = 2Mx 
°dV  @ [A +T2'p)2Sf] 1 Sx (1 +72!) 

When this condition is fulfilled, the transfer function W.gy with respect to the electromagnetic 
torque will completely coincide with expression (2) given in [12]. 

Expression (6) for the link Wpppr has the form of an aperiodic link of the first kind, the gain of which 
in the final form depends on the frequency of the stator supply voltage and on absolute slip. A negative sign 
in front of the fraction indicates that the proposed corrective relationship should be positive. Based on this, 
the designation DPF — dynamic positive feedback was introduced above. In modern drives, it is quite simple 
to realize the necessary dependence of the coefficient on the frequency of the voltage of the electric motor. 

The selection of feedback coefficients according to condition (6) makes it possible to obtain a 
linearized transfer function of an asynchronous electric motor at any frequency of the supply voltage and 
load. Such a transfer function, on the one hand, becomes identical to the previously known formula for the 
transfer function of an asynchronous drive at zero slip (1.e., without load), and, on the other hand, to the well- 
known Kloss formula for static modes of operation of an asynchronous motor. 

The next step shows an example of determining the coefficients of the corrective link at frequencies 
of the supply voltage of 10 and 50 Hz for two types of loads 0.2M,. and M,.. The calculation results of the 
transfer functions of the torque driver of the initial IMD and the corrective link are presented in Table 1. The 
resulting amplitude and phase characteristics are combined with the corresponding characteristics without 
introducing the corrective link (Figure 2a and 2b) and are presented in Figure 4 and Figure 5, respectively. 

In widely used IMDs, it is very difficult to realize the link by mechanical torque. Given that the 
torque is equal to J, Y, and in almost all calculations it is assumed that the rotor flux linkage is constant, 
you can replace the original signal in this local connection with the effective value of the stator current, or its 
active component, which is calculated in all frequency converters. 
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Figure 4. The amplitude and phase characteristics of the IMD at a frequency of the stator supply voltage of 


10 Hz: W1 and W2 - slip at a load moment of 0.2M,. and M,. without correction, W3 and W4 - slip at a load 
moment of 0.2M, and M,.with the introduction of Wppr 


Table 1. The transfer functions of the torque driver of the initial IMD and the corrective element 


id Wimp (P Wppr 
0,038p + 0,226 0,707 

ion Pr 0,0002p2 + 0,0229p + 1,38  0,0038p + 0,226 
á 0,038p + 0,226 3,84 

fz 0,0002p2 + 0,0229p + 4,52 0,0038p + 0,226 
0,027p + 1,548 0,141 

m pi 0,006p2 + 0,628p + 20,56 0,27p + 1,548 
J F 0,027p + 1,548 0,77 
2 


0,006p2 + 0,628p + 21,19 0,27p + 1,548 
* 6, — absolute slip at 0.2M,., 6, — absolute slip at M, 


For stator current feedback, the linearization condition will vary slightly: 
; Ww 
0,8? = —2M,Sy (Tsp + 1) “PPE 


From this equality it follows that the high accuracy of the IMD can be obtained by regulating the 
flow of the rotor, while correcting the linearization conditions. 
However, if there is a permissible design error in determining the parameters of the Wppp link, i.e. 


B  2Mx(T2p+1)-Wy 


sy y, +0 =A 


the amplitude and phase characteristics of the IMD will be similar to the characteristics of a linear link of the 
first kind. 
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Figure 5. The amplitude and phase characteristics of the IMD at a frequency of the stator supply voltage of 
50 Hz: W1 and W2 - slip at a load moment of 0.2M, and M, without correction, W3 and W4 - slip at a load 
moment of 0.2M, and M,.with the introduction of Wppp 


In the last part of the article, the influence of the accuracy of the selection of Wppr link parameters 
on the operation of the IMD is investigated. For this, a transmission coefficient error of + 5% was introduced. 
According to the obtained frequency characteristics, the drive remains stable both at a supply voltage 
frequency of 10 Hz (Figure 6) and at 50 Hz (Figure 7). 
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Figure 6. The amplitude and phase characteristics of the IMD at a frequency of the stator supply voltage of 
10 Hz: W1 - accurate, W2 - when the transmission coefficient of the correction link deviates by 5% 
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Figure 7. The amplitude and phase characteristics of the IMD at a frequency of the stator supply voltage of 


4. 


50 Hz: W1 - accurate, W2 - when the transmission coefficient of the correction link deviates by 5% 


CONCLUSIONS 
In the presented article, the transfer function of the torque former in frequency-controlled IMD s was 


first proposed. This transfer function takes into account the frequency of stator voltage and slip. Such a 
formula, on the one hand, becomes identical to the previously known formula for the transfer function of an 
asynchronous drive at zero slip (1.e., without load), and, on the other hand, to the well-known Kloss formula 
for static modes of operation of an asynchronous motor. 


The frequency characteristics of the IMD, built on this transfer function, gave an adequate explanation 


of the changing dynamics of the drive, noted in experiments conducted earlier. 


The proposed transfer function made it possible to justify the structure of effective correction 


(dynamic positive feedback on the stator current) for sensorless electric drives. 
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